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Abstract 


Over  the  past  three  years  an  extensive  molecular 
laser  modeling  code  has  been  developed  under  ARPA 
sponsorship.  This  model  characterizes  the  time  evolution 
of  pulsed  molecular  laser  systems  including  the  effects 
of  non-Boltzmann  electron  distributions,  stimulated 
emission  processes,  temperature  and  density  transients 
and  state-of-the-art  energy  transfer  rate  coefficients. 
During  the  last  year  a  laboratory  effort  has  been  developed 
to  complement  these  analytic  studies.  This  report  will 
describe  the  results  obtained  from  the  modeling  studies 
and  the  current  status  of  the  laboratory  effort. 


Introduction 


Over  the  past  three  years  an  extensive  molecular  laser 
modeling  code  has  been  developed  at  Wayne  State  University 
under  ARPA  sponsorship.  '  Early  applications  of  this 
dynamic  model  were  made  to  small  signal  gain  measurements 
in  a  TEA  pulsed  CO  laser  system.  The  agreement 

obtained  suggested  the  application  of  the  laser  code  to  new 
CO  laser  systems  in  an  attempt  to  provide  a  guide  to 
experimentation  in  optimizing  the  gain  particularly  on  the 
low  lying  CO  transitions. ^  ^  The  results  of  this  application 
are  in  agreement  with  recent  measurements  in  which  substantial 
9ain  was  found  in  a  low  temperature  pulsed  N-/CO  laser  sys- 
tern. 


Laboratory  experiments  aimed  at  detailed  measurements 
of  small  signal  gain  and  CO  laser  output  are  currently  under¬ 
way  and  an  extension  of  our  laboratory  effort  to  produce  gain 
profiles  and  laser  intensity  versus  wavelength  data  for  use 
in  further  understanding  laser  phenomenology  is  in  progress. 

mhe  quantum  mechanical  close  coupling  code  has  been 
extensively  applied  over  the  past  few  years  to  calculations 
of  the  vibrational  transitional  (VT)  excitation  cross  sections 
for  the  species  C02,H20,H2,C0,N2  and  He.(5)  During  the  last 
year,  this  close  coupling  code  has  been  extended  to  calculate 
vibration-vibration  (W)  cross  sections  for  an  arbitrary 
collision  pair  for  which  general  information  is  available  in 
the  interaction  potential.  Currently  this  code  is  being 
applied  to  CO-CO  W  cross  sections  due  to  the  availability  of 
experimental  data. 

In  the  remaining  sections  of  this  report,  the  current 
results  from  the  molecular  laser  code  are  outlined, 
experimental  data  pertinent  to  our  current  laser  program 
presented  and  a  brief  presentation  of  preliminary  close 
coupling  W  cross  sections  are  given. 
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II.  Molecular  Laser  Modeling  Code 


The  molecular  laser  model  developed  under  the  ARPA/ARO-D 
sponsored  research  program  characterizes  the  time  evolution 
jf  the  separate  vibrational  levels  of  a  multicomponent 
mixture  of  two  molecular  species  and  inert  diluents  in  the 
presence  of  hot  electrons  of  arbitrary  density,  pulse  shape 
and  pulse  duration.  Included  in  this  formulation  are  multi¬ 
quantum  electron-molecule  energy  transfer  processes,  single 
quantum  W  and  VT  processes  and  spontaneous  and  stimulated 
emission  processes.  The  system  of  coupled  rate  equations 
obtained  from  these  processes  has  been  typically  applied 
to  a  system  of  from  20  to  30  vibrational  levels  in  each  of  two 
diatomic  species.  The  time  integration  of  the  resulting 
system  of  coupled  nonlinear  differential  equations  is  performed 
using  a  Runge-Kutta-Merson  procedure  with  an  algebraic  sub¬ 


routine  to  handle  "stiffness".  An  extensive  discussion  of  the 
basic  components  in  this  molecular  laser  model  as  applied 
to  CO  laser  systems  has  been  described  extensively  in  the 
literature  '  and  can  be  referenced  in  Appendices  I  and  u. 

In  the  following  we  outline  the  recent  developments  in  the 
laser  code. 


It  has  been  clear  for  some  time  that  non-Boltzmann  elec¬ 
tron  energy  distribution  functions  generally  exist  in  molecular 
laser  systems.  '  Particularly  as  interest  in  focused  on 
visible  and  UV  excitation  in  molecular  systems,  a  knowledge  of 
the  exact  form  of  the  electron  energy  distribution  becomes 
critical  to  estimates  of  excitation  rates.  In  addition,  the 
interesting  effects  of  gas  additives  on  laser  gain  and  wave¬ 
length  in  CO  suggests  that  the  plasma  properties  are  being  altered 
rather  than  a  direct  effect  via  molecular  energy  transfer 
processes,  (9)  Since  a  shift  in  available  laser  wavelength  toward 
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higher  energy  photons  (lower  vi).  rational  transitions) 

has  been  noted,  particular  interest  is  viewed  toward  plasma 

predictions  in  multicomponent  infrared  laser  systems. 

In  order  to  calculate  the  electron  energy  distribution 
in  the  CO  laser  we  must  solve  the  following  form  of  the 
Boltzmann  equation  for  the  distribution  function  f(u):^ 

2  f  2  2  1 

(1)  E  u  df  +  2m  u  NQ  f  +  kT  u  NQ  df  = 

i  m  m  -  l 

3  NQ  du  M  L  e  du  j 

m  u  ^ 

-nN  ]  f  du'  (u'+u)  0  .  (u'+u)  f  (u'+u) 

^  /  >  3  J 

J  0 

u  . 

-  J  du'  (u' -u)  Q  .  (u^ -u)  f  (u^u)  j 

A  ^  I 


where  u 


m 

M 

N 

E 


T 


electron  energy 

f  h  #  t 

threshold  energy  for  the  j  vibrational  or 

electron  molecular  energy  level 

electron  mass 

molecular  mass 

molecular  density 

electric  fi°ld 

momentum  transfer  cross  section  for  electron- 
molecule  collisions 

cross  section  for  electron  impact  excitation  or 
deexcitation  of  the  jth  molecular  vibrational  or 
electronic  state 
gas  temperature 


The  above  form  of  the  Boltzmann  equation  involves  the 
following  assumptions: 

1)  a  uniform  weak  de-electric  field 

2)  steady  state 

3)  low  density  of  ions,  electrons,  and  excited  molecules 
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4)  large  elastic  cross  section  compared  to  inelastic 
cross  section 

5)  mean  free  path  small  compared  to  the  dimensions  of 
the  laser  tube 

6)  the  proportion  of  ionizing  to  exciting  collisions  is 
negligible . 


If  in  the  form  of  the  Boltzmann  equation  above,  we  make 
the  following  dimensionless  substitutions: 


z  =  eu_  , 
kT 


/ 


g  =  M  |  eE 

6M  \  NQ  kT 


0  =  Q_ 


m 


T)  .  —  M  v  +  j 
2m  ~~ 
Q0 


(where  Qq  is  some  arbitrary  constant  chosen  to  make  the  numbers 
6  and  n  more  manageable,  we  used  Q  =10~16  cm2)  we  obtain 


(2)  Oz^z  +  g/G2)  df  +  0z2f 

dz 


-i  I/’ 


dP(P+z)n.. (p+z) f (P+z) 


/z  . 

Jdr(Z-nr)_.  (Z-P) f (z-P) 


If,  neglecting  the  second  (superelastic)  integral,  we  let 


(3)  g  f  d£(p+z)  .  (p+z)  f  (p+z) 

3 


and  then  differentiate  this  with  respect  to  z  we  obtain 

(4)  d£  =  \  (zfz  )n.  (z+z.)  f  (z+z.)  -  zn.(z)f(z) 
dz^Ll^  3  3  3  J 
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Substituting  (3)  into  (2)  wc  obtain 


( 5  >  df 
dz 


I _  2  zf  +  9  (z)/6z 

z+a/? 


Equations  (4)  and  (5)  are  simultaneous  first  order 
differential  equations  for  the  distribution  function  and  the 
inelastic  collision  integral.  We  solve  these  numerically 
using  a  fourth  order  Runge-Kutta  integration  scheme.  Since 
f(z)  is  strictly  monotonically  decreasing  with  increasing 
energy  z  and  f  (  <»  )  =0  we  start  our  integration  at  large  z  where 
f(z)  is  small  and  integrate  in  toward  zero.  That  is,  the 


distribution  function  is  of  the  form, 
(6)  f  (z)  =  f  K(z')dz" 


00  eo 

= K(z')dz'  -f  K(z')dz' 


■»0  f 

=  JtWT*  +  '  K(z /)  dz 


K(z^) dz  ^ 


Of  course,  instead  of  starting  the  integration  at  infinity 
we  start  at  a  large  but  finite  z  chosen  such  that  the  truncation 
error  is  negligible.  In  this  way  we  are  able  to  deal  with 
the  non-locality  of  the  Boltzmann  equation. 

The  inicial  values  for  f(z)  and  g(z)  are  obtained  using 
an  approximate  solution  to  the  Boltzmann  equation  given  by 
Holstein.  1  This  solution  is  found  assuming  the  following. 
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1)  the  super-elastic  integral  c*n  be  ignored 

2)  the  term  kT  u  NQ  df  which  accounts  for  molecular 

e  du 

motion  can  be  ignored 

3)  the  momentum  transfer  cross  section  is  constant  with 
energy 

4)  the  total  excitation  cross  section  is  a  step  function 
o  for  u  <  u  and  constant  for  u  >  u 

“  G 

a  solution  to  the  Boltzmann  equation  is, 


Ue  =  the  threshold  energy  for  excitation 
Qe  -  the  total  excitation  cross  section 
Qm  =  thG  amentum  transfer  cross  section. 


The  conditions  for  the  validity  of  this  approximation  are, 

2 


usu 


Y»1  ,  2!m 

M 


«  1 


Qp/Qm 

e  m 


Note  that  fH(ue)-l.  That  is,  this  solution  is  unnormalized. 

To  make  this  approximation  a  little  more  realistic  v/e  assume 
that  f(u)  in  the  energy  range  u  £  ue  is  given  by  the  solution 
to  the  Boltzmann  equation  with  the  inelastic  integral  equal 

to  zero.  This  solution  (properly  normalized)  is  Dryv-styn  and 
is 


(8)  f  (u)  =  2^_ 

Hf) 


M  1  ( —\ 

6m  2  \n/ 


m 


.1 

4 


2  2 

u  Qm 
m 


M/3m  (e/n) 


(  i  s  u 


We  now  use  f  (u)  =  fD(ue) *f„ (u) ,  (u>ue>  to  calculate  the 
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initial  values  for  the  Runge-Kutta  integration  sequence. 

This  technique  gives  good  starting  values  but,  as  it  turns 
out,  the  final  result  for  f(z)  after  normalization  is  only 
weakly  dependent  upon  the  starting  values  of  the  integration 
provided  the  integration  is  started  at  large  enough  z. 

Some  typical  results  for  the  electron  energy  distri¬ 
bution  function  using  the  above  approach  are  shown  in 
Figure  1.  The  conditions  for  this  calculation  in  pure  CO  are 
the  same  as  stated  by  Nighan  '  for  comparison.  Nighan’s 
results  are  replotted  in  Figure  la.  Although  qualitatively 
similar,  the  E/N  values  for  the  comparative  calculations  are 
displaced  by  one  curve.  The  source  of  this  discrepancy  is 
under  review.  Clearly,  the  distributions  are  highly  non- 
Boltzmann  and  their  presence  will  strongly  affect  calculations 
involving  excitation  energies  above  about  2  or  3  ev,  depending 
on  the  E/N  value.  The  sensitivity  to  input  cross  section 
values  is  shown  in  Figure  2  where  the  electronic  cross 
sections  are  arbitrarily  reduced  by  a  factor  of  10  (curve  A) 
and  the  vibrational  excitation  cross  sections  reduced  by  2 
(curve  B)  over  the  nominal  values.  The  effect  of  superelastic 
collisions  in  CO  on  the  form  of  the  electron  energy  distribution 
is  shown  in  Figure  3  for  two  CO  vibrational  temperatures 
using  the  nominal  cross  section  values  of  Figure  1. 

The  running  time  for  the  above  calculations  is  approxi¬ 
mately  1  min  CPU  on  an  IBM  360/67  system,  indicating  that  the 
method  is  not  only  stable  but  very  fast.  Since  the  run  time 
for  the  generation  of  the  electron  distribution  function  is  so 
rapid,  it  becomes  feasible  to  perform  a  fully  coupled  electron 
distribution  calculation  with  a  time  evolving  vibrational 
distribution  appropriate  to  a  pulsed  laser  system.  Calculations 
are  currently  underway  to  predict  the  form  of  the  distribution 
function  in  the  case  of  CO  laser  mixtures. 
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Figure  la:  Calculated  electron  distribution  functions  from  the 

publication  of  Nighan(l?)  for  comparison  to  Figure  1. 
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Figure  3:  Variation  in  calculated  electron  energy  distribution 


functions  for  different  CO  vibrational  temperatures 
at  an  E/N  of  5  X  10“  . 


Energy  Transfer  Rate  Coefficients 


Past  applications  of  the  molecular  laser  code  have 
used  W  rate  coefficients  based  upon  the  short  range  force 
model  of  Rapp  and  Englander-Golden^12^  modified  to  fit  the 
experimental  data  of  Hancock  and  Smith. (6,7)  Since  this 
method  is  not  at  all  reliable  for  the  estimate  of  the 
temperature  dependence  of  W  rates  for  near  resonant  processes, 
we  have  added  to  the  laser  code  the  long  range  force  model 
developed  by  Sharmav  ’  and  outlined  by  Jeffers  and  Kelley. (14) 

The  basis  for  this  long  range  force  model  has  been  extensively 

•  •  (14) 

discussed  in  the  literature.  Figure  4  shows  the  present 

W  coefficients  for  CO  at  300°K  as  compared  to  the 

(6) 

available  measurements.  We  are  currently  completing 

the  quadrupole-quadrupole  interaction  case  to  generate  N2~N2 
near  resonant  W  rate  coefficients.  These  will  be  inserted 
into  the  code  for  future  laser  model  predictions . 

Stimulated  Emission  Processes 

Past  published  results  on  the  application  of  the  laser 
code  has  not  explicitelv  included  stimulated  emission  processes 
in  the  prediction  of  the  time  sequence  of  small  signal  gain 
available  from  a  pulsed  CO  laser/  '  '  In  order  to  demon¬ 
strate  this  current  capability,  we  performed  a  calculation 

on  the  room  temperature  TEA  CO  pulsed  laser  used  by 
(3) 

Jeffers  m  which  the  6-5  J=26  (upper  state)  transition 

was  stimulated  and  subsequent  small  signal  gain  calculated. 
Figures  5  and  6  show  the  maximum  gain  profiles  for  a 
particular  vibrational  transition  without  and  with  stimulated 
emission  processes,  respectively.  The  power  output  in  the 
6-5  J=26  transition  is  calculated  to  be  1996  watts  with 
a  pulse  width  (width  at  half  maximum)  of  .23 second. 

A  similar  calculation  for  7-6  J=26  stimulated  is  shown  in 
Figure  7.  The  power  output  for  this  case  is  1047  watts 
with  a  pulse  v;idth  of  ,2  5jX  seconds. 
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As  can  be  clearly  seen,  the  gain  on  both  the  7-6 
and  the  5-4  transitions  are  enhanced  due  to  depopulating 
the  lower  state  and  populating  the  upper  state,  respectively. 
Subsequent  stimulation  of  the  5-4  transition  at  the  optimum 
time  delay  would  produce  enhance  5-4  power  output  over 
that  of  the  unstimulated  case.  The  possibility  of  stimulating 
the  3-2  transition  in  the  cooled  ^/CO  laser  followed 

by  stimulating  the  2-1  transition  could  produce  substantial 
power  output  in  this  high  photon  energy  range.  Calculations 
using  this  stimulated  emission  capability  are  being  performed 
in  an  attempt  to  further  optimize  available  power  output 
on  low  level  transitions. 

Sidelight  Infrared  Spectra 

The  spontaneous  emission  sidelight  spectra  code  is  fully 
operational  either  using  a  predefined  vibrational  distribution 
(i.e.  Treanor  type )  or  the  direct  output  of  the  calculated 
vibrational  distribution  from  the  laser  code.  The  harmonic 
or  first  overtime  spectra  can  be  generated  as  a  line  by  line 
calculation  integrated  over  an  arbitrary  slit  function.  An 
example  of  the  predicted  sidelight  first  overtone  spectra  is 
shown  in  Figure  8  for  a  Treanor  distribution  having  Q*  =2700°K 
and  a  rotational  temperature  of  375°  (the  slit  function  is 
assumed  to  be  triangular  with  a  width  at  half  height  of  6.6  mu . 
A  roughly  equivalent  experimental  spontaneous  emission  spectra 
obtained  from  a  CW  ^/CO  laser  system  in  the  laser  laboratory 
is  shown  in  a  later  section  of  this  proposal. 
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III.  Theoretical  Determination  of  W  and  VT  Rate  Coefficients 

Modeling  of  molecular  gas  laser  systems  requires  a 
knowledge  of  the  reaction  rates  for  a  wide  range  of  excitation 
and  de-excitation  collision  processes  under  a  wide  range  of 
gas  kinetic  conditions.  At  best  this  data  is  difficult  to 
obtain  either  experimentally  or  theoretically,  and  a  combination 
of  approaches  is  necessary  to  estimate  values  for  the  parameters 
involved. 

The  theoretical  program  at  RIES  hac  been  primarily  con¬ 
cerned  with  the  calculation  of  cross  sections  for  the  vibrational 
excitation  of  molecules  by  heavy  body  collisions.  This  is  carried 
out  using  a  numerical  close  coupling  solution  of  the  wave  equa¬ 
tion  for  the  colliding  system  based  on  the  breathing  sphere 
model  of  Schwartz,  Slawsky,  and  Herzfeld  to  represent 

the  target  molecular  vibrator  and  an  empirical  spherically 
symmetric  scattering  potential,  but  neglecting  any  structure 
in  the  impacting  particle. 

To  this  approximation  a  variety  of  collision  systems  have 

*  *  1  *  *  * 
been  studied  (CO  -CO;  CO  -0 {  D) ;  CO,  -  CO,?  CO,  -N, ;  CO,  -  H, ; 

*  *  * 

C02  “  H2°'  H2°  ~  H2°;  H2°  “  H2 '  an<*  H2°  “  C02  w^ere  *  indicates 
the  target  vibrator) .  Many  of  the  results  and  references  are 
collected  in  the  recent  report  by  Fisher  ^  .  However,  these 

studies  were  limited  to  the  case  of  vibrational  -  transla¬ 
tional  (VT)  energy  transfer  processes.  For  many  gas  mixtures 
of  interest  (eg.  CO  -  ^)  direct  vibration  -  vibration  (W) 
ene*-qv  transfer  processes  are  also  important.  Nevertheless 
no  calculations  of  such  processes  to  the  precision  afforded 
by  the  close  coupling  approximation  have  been  published  up  to 
ilo  pr  e  s  e  n  t  t  imo  . 
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The  RIES  molecular  collision  code  has  now  been  extended 
to  include  vibrational  states  in  the  impacting  particle. 

This  in  turn  permits  both  the  calculation  of  significant  W 
cross  sections  simultaneously  with  the  VT  cross  sections 
within  the  same  breathing  sphere  approximation  and  the  estima¬ 
tion  of  the  effect  produced  on  target  particle  VT  cross  sections 
by  coupling  in  the  impact  particle  states. 

Calculations  on  the  CO  -  CO  system  are  currently 
underway  to  debug  the  revised  code,  to  study  the  magnitude 
of  VV  cross  sections  for  states  in  exact  energy  resonance 
and  to  determine  the  effect  of  the  impact  particle  states  on 
the  calculated  VT  cross  section  for  the  excitation  of  the 
first  vibrational  level  of  CO  from  the  ground  state. 

This  cross  section  dominates  the  vibrational  relaxation  of 
CO  and  the  previous  VT  results  were  in  good  agreement  with 
experiment.  Consequently,  for  this  system,  the  VT  cross  section 
should  not  be  particularly  sensitive  to  the  addition  of 
impact  particle  states. 

Preliminary  results  are  shown  in  Table  I  .  Conclusions 
drawn  from  these  values  must  be  speculative  since  the  total 
cross  sections,  which  are  the  physically  important  quantities, 
involve  summation  over  approximately  100  partial  waves  rather 
than  the  single  partial  cross  section  shown.  Nevertheless, 
it  is  interesting  to  note  that  the  partial  cross  section  domin¬ 
ating  the  process  of  vibrational  relaxation  has  been  little 
changed  (30%  decrease)  from  the  original  VT  result  although 
the  VV  cross  sections  themselves  are  large  and  indicate  that 
exchange  processes  suen  aj 

C0(0)  +  CO  (2)  - ►  C0(1)  +  C0(1) 
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are  orders  of  magnitude  faster  in  tne  depopulation  of  tne 
higher  vibrational  levels  than  the  VT  process 

CO  (0 )  +  CO  (2)  —  CO  (  1)  +  (.266  ev)  . 

At  the  present  time  these  CO  calculations  are  being 
extended  over  all  significant  partial  waves,  the  number  of 
coupled  states  required  for  convergence  of  the  cross  sections 
is  being  investigated  and  it  is  intended  to  study  the  sensi¬ 
tivity  of  the  VV  cross  sections  to  the  form  of  the  long 
range  scattering  potential.  Finally,  the  theoretical 
results  will  be  compared  with  the  available  experimental 
values  '  and  a  report  issued  in  the  near  future. 


20 


e 

? 

T 

I 

V 

-C 

£i 

o 

u 

«W 

o 

• 

C-* 
ON  0 
•h  a 
v>  i= 
<0 

4-1  <W 

•H  O 

u 

x  m 

V  4* 

■H 

■h  c 

4  3 

c 

o  c 

•H  -H 
4J  w 

4 

i!  o 

>° 

X  -H 
•U  > 

0  Cl 
•M  > 

«  5 

8* 

•h  ia 

V)  H 
U  4J 
Cl  Vi 

■a 

(0 

V)  Vi 
O  O 
U  TJ 
O  Vi 
0 
•O 

o  o 

4->  VI 
<3  O 
-H  N 


(H  W 

<0  c 
o  o 

■H 

W 

M  -H 

0  W 


c* 

0‘ 

•w 

4J 

o 

« 

M 


4J 

Vl 

«l 


CM  PO 

4> 

— 

, _ , 

00 

00 

H 

O' 

m 

m 

m 

r* 

^  *  1 
H  H  ^ 

1 

1 

i 

i 

1 

i 

i 

i 

i 

*  >91 

CM 

ro 

VO 

VO 

in 

4* 

co 

O  O  CM 

o> 

• 

VO 

• 

H 

00 

O' 

o 

00 

H 

H 

00 

H 

CM 

cn 

CM 

m 

*s 


cm  n 
*  l 


a 

in  to 
«  «' 

Vi 

4J  U 

« 

4  4 

•o 

-V*  4J 

VI 

m  n 

0 

n 

<  a 

V 

u 

0) 

0"0 

N 

C  0) 

H  c 

-1  -H 

a  4 

3  Vi 

O  0) 

Li  U 

m  -*  ^ 

00 

•H  ^  O'  P'1 

1 

w  «L  >L  J, 

CM 

p^  m  vo  cm 

vo  oi  »*. 

•  , 

•  •  •  • 

«e 

H  n  (M  n 

00 

I 


VO 

VO 


(0 

m 

« 

8 

a4 

c 

o 

■ri 

in 


»* 

>  c*> 
Cl  p* 
H 

♦  -f 
— —  — 
O  ^ 

C  o 

V  V 


>  > 
<D  Cl 
00 

<n  to 


5 

»• 

cn 


> 

O 

oo 

VO 


>  > 
0)  01 
00  M* 
VO  PO 


♦  + 

■ — -  .  ^ 

o  rvi 

O  O 

V  V 


+  ♦ 

rH  CM 

o  o 

V  V 


+  + 
fM  O 

o  o 

V  V 


®  ®  '  I  ®  ^  H  ri  n  u 

o  o  ~  ~  ~ - 

ci  ci 


> 

o 


> 

Cl 


©  O 

o  o 

V  ? 

o  c 


Cl  Cl 

« 


o 

rH 

rH 

CM 

»H 

•H 

O 

n-* 

>— *■ 

w 

•-.r- 

o 

o 

O 

o 

o 

o 

O 

u 

0 

o 

Cl 

Cl 

o 

C) 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

♦ 

> 

> 

> 

> 

> 

> 

Cl 

Cl 

Cl 

Cl 

o 

o 

00 

4f 

00 

00 

> 

ro 

ro 

Mi* 

rn 

VO 

V0 

Cl 

H* 

4J* 

• 

• 

• 

• 

• 

JL 

X 

+ 

+ 

X 

4 

_ _ _ 

o 

•H 

•H 

CN 

rH 

CM 

CN 

w 

o 

o 

O 

o 

o 

o 

o 

V 

y 

y 

u 

y 

u 

u 

4- 

4 

4 

+ 

4 

+ 

p— * 

o 

rH 

rH 

CN 

o 

O 

O 

■V^ 

-M* 

W 

- 

- 

O 

0 

0 

u 

Cl 

u 

o 

§ 

4J  CM 


h  gim 

a.c  o 

■H  ±1 
•V  4J 

H  VO 
3H  ti 

e  <m 

C  «M 
V  0)  0) 

Ja  «o 

■ri  0> 

o  >  j: 

*»  w  +i 

0)  •  10 
3  0  0) 
•H  U  -P 

*  4 

O 

..  OH 

sz  'o 

o  c  c 

•«H  O-H 
■C  -H 
>  4J  <U 

4  o 

>1  X  c 

•Q  4  Cl 
•H  Vi 
cl  m  in 
01  Vi  vu 
4J  IM 


3 
eo 
01 
Vi 

I 

> 

to 
3 
•  0 
0)  -H 
> 
*4  01 
0  Vi 

<w  a 

O'  0) 
C  £ 
•H  4J 
•H 

CUJZ 
3  +1 
0-H 
f»  » 


>  > 
0)  0) 
00  4T 

VO  rr, 
M1  f-> 


V 


<M  4  Q 
O  C 

o  • 

V-HH 
0)  4-1  10 
3  4  U 

a  Vi  -H 
-Q  4-< 
-C  C 
«  >  0) 
0)  13 

VI  01  H 
4  .C 

3*1! 

TJ  « 

C  0)  -o 

•H  41  H 
4  J 
«  C  O 
h-hj; 
«  e  in 

.X  0 
y  T3  in 
4  C 
Vi  to  o 
.O  in  -ei 

<U  4J 

s  y  o 

••H  w  V 

Vi  M 

©  a 

Vi-  « 

3  in 

OV  H  o 

•H  JC  14 

fci  c-<  y 


0)  r-i 

4 
4  U 

V>  -H 

n  4J 

c 

i*>  0) 

•o 

+>  -H  • 

in  -a  • 
c  >■  a;  in 

•h  —i  u  o 
4  vi  u  in 
D<  ill  o  in 
4  0(PH  Cl 

o  O'  a 

W  VI  a*  O 

4  a  c  vi 

a 

<  «  < 

e 

n  o  c 

O  p 
<M  4-1 

u  o  <u 

O'  O  Vi 

c  in  3 

•*4  4-1 

-c  W  u 

a  in  3 

3  0 


> 

I 

I  > 


o 

c 

4 

c 

o 


o  o  m 


Q)  Cl  ••  4J 

vi  m  .  U 
4  Cl  Cl  4 

4J  x:  .  x 

in  -i  y 


21 


IV.  Laser  Experiments 


In  this  section  we  outline  the  progress  made  in  the 
laser  experiments  over  the  past  year.  Our  effort  has  been 
aimed  at  providing  details  of  the  gain  profiles  and  vibra¬ 
tional  distributions  present  in  CW  and  pulsed  CO  laser 
systems  for  use  in  conjunction  with  modeling  studies  to 

completely  elucidate  the  phenomenology  operative  in  these 
laser  systems. 

The  master  laser  oscillator  was  put  into  operation 
during  the  summer  of  1972.  Basically  the  system  was  composed 
of  an  8  nun  diameter  discharge  water  cooled  tube.  The 
discharge  was  run  from  the  ends  of  the  tube  to  the  center. 

An  enhanced  metal  coated  stainless  steel  mirror  was 
used  at  one  end  and  a  grating  wavelength  selector  in  Littrow 
configuration  was  used  as  the  reflector  at  the  other  end. 
Difficulties  were  encountered  in  room  temperature  operation. 

A  variety  of  conditions  was  tried  copying  most  of  the  published 
data  but  no  output  was  measured.  The  tube  was  then  cooled 
to  dry  ice  temperature  and  CW  lasing  was  monitored  on  more 
tnan  forty  CO  spectral  transitions.  Powers  ranged 
from  microwatts  to  tens  of  milliwatts. 

A  small  amount  of  ripple  Is  present  in  our  DC  power 
supply.  This  results  in  large  (up  to  50%)  intensity  fluctu¬ 
ations  in  the  laser  output.  An  attempt  made  to  filter  this 
using  an  RC  filter  configuration  was  partially  successful  but 
it  was  felt  higher  stability  was  required.  A  current  regu¬ 
lator  was  constructed  after  the  design  of  Posakony  (17>  and 
this  limited  the  output  fluctuations  to  no  more  than  a  few 
percent.  The  probe  oscillator  was  then  ready  for  operation. 

The  test  discharge  cavity  is  not  year  operational  due  to 
difficulties  with  the  suppliers  of  the  gas  monitoring 
equipment.  In  the  meantime  while  awaiting  supply,  it  was 
decided  to  attempt  to  understand  why  room  temperature  operation 


\ 
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could  not  be  achieved.  Conditions  in  the  discharge  tube 

were  monitored  via  end  view  spontaneous  emission  measurements 

under  a  variety  of  operating  conditions.  Analysis  of  this 

data  is  proceeding  from  two  directions.  in  one  effort 

artificial  spectra  are  being  generated  in  an  attempt  to 

match  the  observed  spectra.  in  the  other  method  a  least 

squares  fitting  program  has  been  written  which  combines 

gradient  search  techniques  with  an  analytic  solution 

developed  from  linearizing  the  fitting  function.  This 

method  has  been  outlined  by  Bevington ^ ,  and  follows 

(19) 

the  technique  of  Marquardt  .  The  program  is  operational 
and  when  tested  on  an  artificially  generated  spectrum  it 
obtained  the  vibrational  temperature  to  within  0.7%  in  a 
time  of  about  3  seconds. 

Figure  9  displays  the  intensity  versus  wavelength 
out Dut  from  the  grating  selected  master  oscillator  system 
at  dry  ice  temperature.  As  noted  by  Bhaumik,  et.  al.^20\ 
upon  the  addition  of  Hg  vapor,  the  spectral  output  shifts 
as  shown  in  Figure  10.  The  details  associated  with  this 
shift,  and  possible  shifts  induced  by  the  presence  of 
other  gas  additives,  may  be  very  important  in  optimizing 
CO  laser  output.  Experimental  as  well  as  analytical 
efforts  are  currently  aimed  at  understanding  this  enhanced 
laser  output . 

A  typical  spontaneous  emission  endview  measurement  is 
shown  in  Figure  11  for  dry  ice  conditions.  Preliminary 
interpretation  of  the  spectra  shows  that  the  rotational 
temperature  is  considerably  higher  than  the  wall  temperature 
which  undoubtedly  is  the  reason  that  no  lasing  is  observed 
under  room  temperature  operation. 
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